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ABSTRACT 

The high densities, long hfetimes, and narrow emission measure distributions 
observed in coronal loops with apex temperatures near 1 MK are difficult to rec- 
oncile with physical models of the solar atmosphere. It has been proposed that 
the observed loops are actually composed of sub-resolution "threads" that have 
been heated impulsively and are cooling. We apply this heating scenario to nearly 
simultaneous observations of an evolving post-flare loop arcade observed with the 
EVYl/ STEREO, XRT/Hmode, and TRACE imagers and the EIS spectrometer 
on HINODE. We flnd that it is possible to reproduce the extended loop lifetime, 
high electron density, and the narrow diflFerential emission measure with a multi- 
thread hydrodynamic model provided that the time scale for the energy release 
is sufliciently short. The model, however, does not reproduce the evolution of 
the very high temperature emission observed with XRT. In XRT the emission 
appears diffuse and it may be that this discrepancy is simply due to the diffi- 
culty of isolating individual loops at these temperatures. This discrepancy may 
also reflect fundamental problems with our understanding of post-reconnection 
dynamics during the conductive cooling phase of loop evolution. 

Subject headings: Sun: corona 

1. Introduction 

One of the principal problems in solar physics is understanding how the Sun's corona 
is heated to very high temperatures. Recent work on coronal loops indicates that they have 
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physical properties that are difficult to reconcile with theoretical models. Coronal loops 
with temperatures near 1 MK are observed to persist longer than a characteristic cooling 



time, suggesting steady or quasi-steady heating (e.g., iLenz et al.l Il999l : lAschwanden et al 



2000). Steady heatin g models, however, can not reproduce the high electron densities ob- 
served in these loops (jWinebarger et al.ll2003l ). Mul ti-thread, impulsive heating models have 



been proposed as a possible heating scenario (e.g., ICargill fc KlimchukI 119971 : 1 Warren et al. 



20031 : IPatsourakos fc KlimchukI l2008l ) . Such models are motivate d by our und erstanding 
of the energy release during magnetic reconnection in flares (e.g., IParkerl 1 19831 ). In these 
models impulsive heating leads to high densities and multiple, sub-resolution "threads" 
lead to long lifetimes relative to the cooling time for an individual loop. These models 
are severely constrained by the relatively narrow distributions o f temperatures that are of- 
ten observed in loops with apex temperatures near 1 MK fe.g.. IPel Zanna fc Mason 2003 : 
Aschwanden fc Nightingale! [20051; ICirtain et al.l l2007l : I Warren et al.l l2008l ) . A narrow distri- 
bution of temperatures suggests that the loop can contain only a few independent threads. 

One difliculty with fully testing coronal heating scenarios such as these with hydrody- 
namic models has been the spareness of data. Previ ous work on loop evolution has generally 
focused on measurements imaging instruments (e.g.. Warren et a02003 : Aschwanden fc Nightingale 



20051 : lUgarte-Urra et al.ll2006l ). which have limited diagnostic capabilities. Current solar ob- 
servatories, however, allow for coronal loops to be observed in unprecedented detail. The 
EUV Imaging Spectrometer (EIS) on the Hinode mission provides high spatial and spectral 
resolution observations over a very wide range of coronal temperatures. EIS plasma diag- 
nostics yield important constraints on the physical properties of coronal loops. The X-ray 
Telescope (XRT) on Hinode complements these observations with high spatial and tempo- 
ral resolution observations of the high temperature corona. The multiple viewpoints of the 
twin STEREO spacecraft allow for loop geometry, a critical parameter in the modeling, to 
be measured using the EUV Imagers (EUVI). The Transition Region and Coronal Explorer 
(TRACE) currently provides the highest spatial resolution images of the solar corona. 

In this paper we use STEREO^ Hinode^ and TRACE observations of an evolving loop in a 
post-flare loop arcade to make quantitative comparisons between a multi-thread, impulsive 
heating model and measured densities, temperatures, intensities and loop lifetimes. An 
important component of this work is the development of methods for integrating the different 
observations into hydrodynamic simulations of the loop. We flnd that it is possible to 
reproduce the extended loop lifetime, the high electron density, and the narrow differential 
emission measure (DEM) with a multi-thread model provided the time scale for the energy 
release is sufficiently short. The model, however, does not reproduce the evolution of the 
high temperature emission observed with XRT. 
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One goal of investigating the heating on in dividual loops is to motivate the modehnR of 
entir e active regions o r even the full Sun (e.g ., ISchrijver et al.ll2004l : IWarren fc Winebarger 



20061 : iMok et al.l l2005l : iLundquist et al.ll2008l ). It is possible, however, that there is not a 



single coronal heating mechanism that can be applied to all coronal loops. For example, it 
may be that steady heating is the dorainant heating scenario on some fraction of coronal loops 
(e.g., lMartensll2008l : lAntiochos et al.ll2003l ). Even if impulsive heating of the kind discussed 
here is only a minor contributor to the heating of the solar corona, this study provides 
important insights into the energy release during magnetic reconnection, a fundamental 
process in astrophysical and laboratory plasmas. 



2. Observations 

In this section we provide an overview of the instruments and observations used in this 
study. A summary of the observations is shown in Figure [H The loop considered here is a 
post-flare loop from a very small event (GOES class B2.5) that peaked around 19:00 UT on 
May 2, 2007. 

The EIS instrument on Hinode produces stigmatic spectra in two wavelength ranges 
(171^212 A and 245-291 A) with a spectral resolution of 0.0223 A. There are V' and 2'' slits 
as well as 40'^ and 266'^ slots available. The slit-slot mechanism is 1024^' long but a maximum 
of 512 pixels on the CCD can be read out at one time. Solar images can be made using one 
of the slots or by stepping one of the slits over a region of the Sun. Telemetry cons t raints 



gen erally limit the spatial an d spectral coverage of an observation. See lCulhane et al.l ( 120071 ) 



and iKorendyke et al.l (120061 ) for more details on the EIS instrument. 



For these observations the slit was stepped over the active region and 15 s exposures 
were taken at each position. An area of 256^' x 256^' was imaged in about 71 minutes. A total 
of 20 spectral windows were read out of the CCD and included in the telemetry stream. The 
raw data were processed using eis_prep to remove dark current, warm pixels, and other 
instrumental eflFects using standard software. During the processing the observed count 
rates are converted to physical units. Intensities from the processed data are computed by 
fltting the observed line pr oflles with Gaussia ns. The EIS rasters are co-aligned to account 



for any spatial offsets (see lYoung et al.ll2009l for a discussion). Spacecraft jitter during the 



raster has not been accounted for. The Hinode housekeeping logs suggest relatively small 
displacements (less than one pixel) for the narrow fleld of view of interest here. For larger 
structures spacecraft jitter can be important. EIS rasters in a number of diflFerent emission 
lines are shown in Figure O and show post-flare loops at various temperatures in the lower 
part of the active region. These rasters also indicate a brief data gap due to orbital eclipse. 
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One limitation of these EIS data is the lack of temporal information. Better information 
on the tem poral evolution of t hese loops is provided by the imaging instruments, such as 
the EUVI faoward et aPboOsI ) on the Solar Terrestrial Relations Observatory (STEREO) 
mission. The EUVI is a normal incidence, multilayer telescope which can observe the Sun 
in 4 wavelength bands centered at 284, 195, 171, and 304 A. EUVI observes the full Sun and 
therefore has reduced spatial resolution (1.6^^ pixels) relative to the other observations that 
we consider here. There are two STEREO spacecraft with identical instrument packages. 
The twin STEREO spacecraft drift away from the Earth at about 23° per year. On May 2, 
2007 the separation between the spacecraft was small, about 6°. 

The EUVI images taken around the time of the EIS raster are indicated in Figure [H 
Because of telemetry constraints the image cadence is limited. For these observations 171 A 
images were taken at a relatively high cadence 150 s) while the images at the other 
wavelengths were taken at lower cadences 600-1200 s). The raw data are processed using 
euvi_prep to produce calibrated, co-aligned images. Images of the active region and flare 
from STEREO B EUVI are shown in Figure [31 

The XRT on Hinode is a high cadence, high spatial resolution (approximately V' pixels) 
grazing incidence telescope that images the Sun in the soft X-ray and extreme ultraviolet 
wavelength ranges. Temperature discrimination is achieved through the use of focal plane 
fllters. Because XRT can observe the Sun at short wavelengths, XRT images can observe high 
temperature solar plasma very efliciently. The thinner XRT fllters allow longer wavelength 
EUV emission to be images and extend t he XR T response to lower temperatures. Further 



details on XRT are given in iGolub et al.l (120071 ). 



As indicated in Figure [U the principal XRT images taken during this time period were 
in three fllters, Ti-Poly, Al-Thick, and Be-Thick, at a variable cadence. Unfortunately, the 
exposure times on the Be-Thick images are too short for the images to be used for analyzing 
active region loops. The standard processing routine xrt_prep is used to remove the CCD 
bias, dark current, and calibrate the images. Images are also "dejittered" using Hinode 
spacecraft housekeeping data so that the images are co-aligned with respect to each other. 
Hinode tracks solar rotation so there is no need to account for it in the images. Example 
Ti-Poly images are shown in Figured! These images have a fleld of view of 512^' x 512'^ 

The TRACE instrument is a high resolution normal incidence telescope. The primary 
and secondary mirrors are divided into quadrants and a rotating shutter is used to select 
which quadrant is illuminated. Three of the quadrants are coated with multilayers for imag- 
ing at EUV wavelengths. The multilayer coatings have peak sensitivities at approximately 
171, 195, and 284 A. The fourth quadrant is coated with aluminum and magnesium fluoride 
for imaging very broad wavelength ranges near 1216, 1550, 1600, and 1700 A. Ima eces in all 



-5- 



of the wavelengths are projected onto a single detector, a 1024 x 1024 CCD. Each CCD pixel 
re presents a solar area approximately 0.5^' on a side. The instrument is described in detail 
by iHandv et al.l (119991). T he initial in-flight performance is reviewed by IColub et al.l (119991 ) 
and ISchrijver et al.l (119991 ) . 

During this period TRACE observed mainly in the 171 A channel at a cadence of about 
60s with occasional 1600 A and white light context images. All of the images have a 512^'' x 
512^' fleld of view. As shown in Figured! there are periodic data gaps in the TRACE data 
due to orbital eclipses. All of the TRACE images are processed using a standard application 
of trace_prep. Additionally, the images are despiked and co-aligned with respect to each 
other to account for solar rotation and drifts in the pointing. A simple cross-correlation 
method is used for this purpose. The TRACE 171 A images are very similar in appearance 
to the EUVI 171 A images. 



3. Analysis 

The primary goal of this study is to compare multi-thread, hydrodynamic simulations 
with the emission observed in an evolving coronal loop. As we will discuss in more detail later 
in the paper, hydrodynamic simulations involve solving the equations for the conservation 
of mass, momentum, and energy in the loop given some input heating rate. Relating the 
heating rate to physical observables is a critical element of the modeling. Our strategy is 
to measure the electron density with EIS and use a family of hydrodynamic simulations to 
infer the required heating rate for this density. Previous numerical simulations suggest that 
for a fixed loop length there is a power law relatio nship between the peak electron density 



and the input energy (IWarren fc AntiochosI l2004l ) . Since we want the loop length to be 



fixed, the other critical element of this modeling is an accurate measurement of the loop 
geometry, including the inclination. Observations from the twin EUVI instruments allow 
the loop geo metry to be measured, and we use the STEREO software package developed for 



this purpose (lAschwanden et al.ll2008al ). 



Once the density and loop geometry have been determined we can perform hydro- 
dynamic simulations and synthesize the expected emission during the entire evolution of 
the loop. The simulation results can then be compared with light curves determined from 
TRACE and XRT and in this section we discuss how these light curves are calculated. The 
distribution of temperatures in the loops is an important constraint on the modelings and in 
this section we also discuss the calculation of the differential emission measure distribution 
with the EIS spectra. 
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3.1. EIS 

Our analysis requires the identification of loops observed simultaneously with both EIS 
and TRACE. To facilitate this we wrote routines to co- align EIS rasters with TRACE images 
and to display 24 bit color images using an EIS raster for one color channel and a TRACE 
image for another color channel. An animation of these images allowed us to quickly identify 
times when the EIS slit was co-spatial with a loop observed with TRACE. TRACERS small 
field of view and frequent data gaps due to orbital eclipse make finding good data sets more 
difficult than anticipated. 

To optimize the co- alignment between EIS and TRACE we had to allow for a roll angle 
between the images. This is in addition to the usual spatial shifts between the pointing 
information contained in the data headers. Since the TRACE data was taken in the Fe IX/x 
171 A channel we used the EIS Fe X 184.536 A raster for co-alignment. 

The EIS intensities for the loop of interest are summarized in Figure [51 Once the region 
of interest was identified we manually selected spatial positions along the loop in the EIS 
Fe XII 195.119 A raster. These points are used as spline knots to define the loop coordinate 



system with s along the loop and t perpendicular to it (see lAschwanden et al.ll2008a 

Figure 3). Since these loop coordinates are not necessarily aligned to the CCD we have 
interpolated to determine the intensities along the selected segment. The loop segments 
displayed in Figure [5] are interpolated to 0.2^' per pixel. Using the loop coordinate system it 
is a simple matter to compute the intensity averaged along the loop segment. The coordinate 
system derived from Fe XII 195.119 A is used for all of the EIS rasters. 

To further isolate the intensity of the loop we identify two background points and fit 
a single Gaussian with a linear background to the selected region. Background subtraction 
is essential to separating the intensity in the loop from the contribution of the ambient 
corona, but there is no unique method for computing it. Analysis of EUVI data, which 
has the advantage of providing two different lines of sight for a single loop, suggests that 
the background subtracted intensities can be computed consistently, although the re can be 



considerable uncertainties for individual measurements (lAschwanden et al.ll2008al ) 



Since we are interested in emission that is co-spatial we extract the same region from 
all of the EIS rasters. The resulting intensities are shown in Figure [51 To test how co- 
spatial the emission is at various temperatures, we have calculated a simple correlation 
coefficient between the background subtracted intensity in each line and Fe XII 195.119 A, 
which represents a middle ground between the highest and lowest temperature emission that 
is observed. For this loop well correlated, co-spatial emission is observed for Fe XIII and 
below. For the emission at higher temperatures (Fe XIV-xvi) the correlation is poor. 
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It is clear from the images shown in Figures [2] and [5] that the combination of high spatial 
resolution and good temperature discrimination allows EIS to probe the interrelationship of 
emission at different temperatures. The loop intensities suggest that the emission at different 
temperatures is generally not co-spatial and that the DEM in the loop of interest should 
be relatively narrow. This is shown more clearly in the "multicolor" image of this region 
presented in Figure [61 This image, which is a 24 bit image formed from rasters in 3 different 
emission lines, would be white in regions where the DEM is broad and the emission is strong 
in all three lines. There are some composite regions that are cyan (green + blue), but the 
post-flare loop arcade is generally dominated by emission in the primary colors suggesting 
relatively narrow distributions of temperature in each loop. 

To investigate the temperature structure of this loop more quantitatively we compute 
the differential emission measure using the background subtracted loop intensities for the 
loop segment. The intensities are related to the differential emission measure by the usual 
expression 

h = ^ [ ex{n,^T)aT)dT, (1) 



47r „ 

where ex{ne^T) is the plasma emissivity and ^{T) is the differential emission measure distri- 
bution. We consider a Gaussian DEM model 



U^) = 7^ exp 



{T - Tof 
2al 



(2) 



which allows for a dispersion in the temperature distribution. Since the density is an im- 
portant parameter in determining the emissivities of many of these lines we leave it is a free 
parameter. To determine the best-fit parameters {EMq^ Tq, a^, n^) we use a Levenberg- 
Marquardt tec hnique implemente d in the MPFIT package. The CHIANTI 5.2 atomic physics 



database (e.g.. Landi et al. 20061 ) is used to calculate the emissi vities. The abundance s of 
Feldman et al. and the low density ionization fractions of iMazzotta et al.l (119981 ) are 



assumed. 

There are several subtleties to computing the emission measure parameters. One is the 
statistical uncertainty associated with each intensity. Because the intensities are averaged 
over a significant area the statistical uncertainties are generally small. The systematic errors 
introduced by the background subtraction and the atomic data, however, are large, but 
difficult to estimate. We simply assume that the relative error is 20% of the measured 
intensity. Another question is how to deal with the emission from lines such as Fe XVI 
262.984 A, that do not show evidence for the loop and have no measured intensity. For 
these lines the background subtracted intensity is zero. Observations of these lines provide 
important constraints on the high temperature component of the DEM and must be included. 
To account for possible errors in computing the background subtracted intensities for these 
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lines we assume that the uncertainty in the intensity is 20% of the background instead of 
20% of the loop intensity. Since the background can be large this represents a substantial 
enhancement of the uncertainty. 

The resulting DEM for this loop segment is shown in Figure [51 The observed and 
computed intensities are given in Table [H From this analysis we obtain an electron den- 
sity of logne = 9.7. For this work we use the Fe XIII 203.826 /202.044 A lines to provide 
the bulk of the density sensitivity. Recently Young et al. ( 20091 ) have noted systematic dis- 



crepancies between the various density sensitive line ratios from Fe XII and Fe XIII. In 
light of this we compared the densities inferred from Fe XIII 203.826/202.044 A and Fe XII 
186.880/195.119 A with those from Si X 258.375/261.058 A in a series of other active region 
and quiet Sun observations. The Si X lines, which were not included in this study, are rela- 
tively weak and not sensitive over as large a range as the Fe XII and Fe XIII lines. However, 
the atomic data for Si X is potentially more reliable than the atomic data available for the 
complex Fe ions. We find that the densities derived from the Si X and Fe XIII ratios are 
in excellent agreement and emphasize the Fe XIII ratio here. The densities inferred from 
FeXII 186.880/195.119 A are as much as a factor of 3 higher. Details of this analysis will be 
presented in a future paper. 

The dispersion in the temperature is found to be logar = 5.5, w hich is comparable 



to other active region loop observations with EIS ( j Warren et al.l l2008l ) . These lines ob- 



served in the quiet corona above the l imb indicate much narrower temperature distributions 



(logcTy ^ 5.0. IWarren fc Brooksll2009l ). Here we find a dispersion in temperature that is sev- 
eral times greater, indicating that this loop is not strictly isothermal. The relatively intense, 
co-spatial emission observed from both Si VII and Fe XIII provides the best direct evidence 
for a distribution of temperatures in this loop. These lines have peak temperatures of for- 
mation that are about IMK apart. The application of a delta function emission measure 
to Equation [1] confirms that a single temperature model cannot adequately reproduce the 
intensities observed in these emission lines. 



3.2. TRACE 

In searching for coronal loops observed with both EIS and TRACE we co- aligned the EIS 
rasters with the TRACE images for this period. Thus to compute the TRACE intensities 
we use the same coordinate system and apply it to all of the TRACE images taken during 
the time of interest. Since the co- alignment between EIS and TRACE is not perfect, the 
spline knots selected in the EIS raster are modified slightly to better align with the loop. 
Background subtracted intensities are computed for each of the available TRACE images 
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taken during this time using the same procedure that was used on the EIS data. For each 
image the region defined by the sphne knots was extracted, straightened, and averaged along 
the loop coordinate to produce the intensity as a function of the perpendicular coordinate. 
The intensity at each time refers to the background subtracted intensity integrated over the 
loop. 

The TRACE light curve is shown in Figure [71 For most of this period the background 
subtracted loop intensities are in the noise and the evolution of this loop can be seen clearly. 
We have manually selected the region around the peak in the light curve and fit it with a 
single Gaussian. This fitting yields a Gaussian width of about 294 s. 



3.3. EUVI 

To extract the three dimensional geometry of this loop we use the software package 
developed by Markus Aschwanden and t he STEREO team f or this purp ose. An application 



of this software is discussed in detail in I Aschwanden et al.l (j2008bl ) and lAschwanden et al 



(l2008al ). The initial processing co-registers images from STEREO A and B to account for 
differences in spacecraft roll angle and spatial resolution. The next step is to outline the loop 
in the STEREO A image. The selected coordinates are projected onto the corresponding 
B image. Since geometry of the loop has not yet been determined this projection yields a 
range of possible coordinates in the B image and the user selects the position of the loop 
within this range. Once the coordinates of the loop have been selected in both images the 
geometry of the loop is determined from simple trigonometric relationships. 

Due to the gap in the A data, an image pair at 21:16 UT, after the peak emission 
observed in this loop, is used, so there is some ambiguity in the loop identification. The 
selection of nearby structures in the loop arcade generally yield very similar results. The 
projection of the extracted loop geometry onto the EUVI B image available at the peak of 
the loop emission (20:59 UT) also outlines the observed loop very well. Thus it does not 
appear that the calculation of the loop geometry is significantly impacted by the data gap. 
The selected loop and the projection of this loop in various planes is shown in Figure [51 



3.4. XRT 

Images taken with the XRT provide information on the evolution of the loop at high 
temperatures. We have computed an XRT light curve similar to that computed for TRACE. 
To use the loop coordinates derived from EIS we first co-align the XRT images with the 



EIS Fe XVI 262.984 A raster, which shares some common features. We then compute the 
background subtracted intensities by selecting two background points and doing a hnear fit. 
The emission seen earlier in the event is much broader than what is observed TRACE and 
so we select a wider area to compute the intensities. The resulting XRT light curve is shown 
in Figure M 

As is indicated by the light curve, the XRT images clearly show strong emission in 
the region that is eventually occupied by the loops observed with EIS and TRACE. It is 
also clear, however, that XRT does not show any individual loops that are as narrow as 
those that are seen at cooler temperatures. Consistent with this, the loop cross section 
measured with XRT is systematically wider than what is measured with EIS and TRACE. 
To illustrate these differences we have constructed multicolor images from XRT and TRACE. 
These images, which are presented in Figure [101 ^se different color channels to display the 
XRT and TRACE data in the same picture. To illustrate the differences in morphology as the 
plasma cools we have offset the times of the selected images by one hour. We show examples 
of the loop cross sections in Figure [TTl This difference between high tempe rature and low 



temperature emission in flares and active re gions is well documented (e.g., iTripathi et al 



20091 : IPatsourakos et al.l l2002l : IWarrenI l2000l ) , and may be further evidence for flne scale 



structure in the solar corona. We will discuss this point in more detail later in the paper. 



4. Hydrodynamic Modeling 

One of the primary paradoxes of coronal loops with temperatures near 1 MK is the 
disparity between the rapid cooling suggested by the high electron densities and the relatively 
long observed lifetimes. EIS density diagnostics allow us to make more rigorous comparisons 
between these time scales. Densities inferred from observed intensities requires an accurate 
determination of the differential emission measure, the absolute instrumental calibration, 
and the loop geometry. This measurement also represents a lower bound on the density. 
The density inferred from density sensitive line ratios circumvents many of these problems. 

If we assume that the loop is cooling only through radiation the energy balance is simply 

— ^ -nlA{T^), (3) 

where A(Te) is the radiative loss function for an optically thin plasma. In the limit of no 
flows the energy is = |P. The plasma pressure given by P = 2ne/c^Te, with ks is the 
Boltzmann constant. The radiative cooling time is deflned as 
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and is given by 



neA{Te) 

Using the temperature and density derived from the EIS DEM a nalysis (loRn^ 9 . 7 and 
logTg :^ 6.11) and a radiative loss rate of 3.2 x 10~^^erg cm^ s~^ ( IBrooks fc WarrenI l2006l ) 
we obtain a radiative cooling time of 341 s. 



The radiative cooling time is not directly comparable to the loop lifetime that we have 
measured with TRACE. If we make the additional assumptions that Te{t) = e'x.\){—t / tt) 
and Te(t) ^ Ti^{t) we can relate the radiative cooling time {tr) to the timescale for changes 
in the temperature (r^) 

Tt = Tr. (6) 

a 



Numerical simulations suggest that a ^ 2 (jjakimiec et al.lll992l ). Finally, to compare with 
the observed loop lifetime we must incorporate the temperature changes into the TRACE 
temperature response curve. This yields 



Ix{t) - exp 



2t2 



(7) 



where T\ is the peak tem perature of the TRA CE response and ax is the Gaussian width of 
the TRACE response (see lWarren et al.ll2003l Equation 9). This yields 



a\l + a 



-Tr. 



With TiTi = 0.96 MK and am = 0.25 MK flWarren et al.ll2003h we obtain at 
is smaUer than the observed Gaussian width of 294 s. 



(8) 



141 s, which 



This mismatch between the predicted and observed loop lifetime is one of the key mo- 
tivations for the multi-thread modeling of coronal loops. By assuming that the observed 
emission comes from a series of loops that are heated at different times it is clear that we 
can create a composite loop with the required lifetime. The challenge is to also match the 
relatively narrow DEM and the lifetime of the loop as it is observed with XRT. 

To simulate the evolution of this loop we consider numerical solutions to the full hy- 
drodynamic loop equations using the NRL Solar Flux Tube Model (SOLFTM). We adopt 
many of the same parameters and assumptions that were used in previous simulations with 
this code an d we refer the reader to the earlier p apers for additional details on the numerical 



model (e.g., iMariskal 119871 : iMariska et al.lll989l ). For example, we assume that the loop is 



symmetric and only simulate the evolution over half of the loop length. We also assume that 
each loop has a constant cross section. 
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For this work we consider a heating function for each thread that is a simple spatiaUy 
uniform heating rate 



where g{t) is a triangular envelope that peaks at time U and has width 5^ which we set to 
be 100 s. The parameter Eq is a small background heating rate that provides the initial 
loop equilibrium. As we mentioned earlier, our strategy for inferring the heating rates from 
the observations is to use the density determined from EIS and the results from systematic 
hydrodynamic simulations. This will determine the peak heating rate for the ensemble of 
threads. The heating rate for the other threads will be determined by assuming a Gaussian 
envelope for the heating function. 

To determine the relationship between the input heating rate and the densities and 
temperatures observed during the cooling phase of the loop evolution, we have performed 
21 simulations with varying between 10~^ and 10^ erg cm~^ s~^. Using the parameters 
derived from the EUVI observations, the loop length is fixed at 135 Mm and the loop in- 
clination is fixed at 68.5°. For each simulation we average over the loop apex to determine 
a representative density and temperature. The simulation results are summarized in Fig- 
ure [121 We find that the apex density at 1.3 MK, which is the peak temperature in the 
DEM, essentially scales as ^ V^f- This implies that for a fixed loop length the observed 
intensity is linearly proportional to the input energy. The relationship shown in Figure [l2] 
indicates that a heating rate of 0.8 erg cm~^ s~^ is required to reproduce an apex density of 
log rig c:^ 9.7. 

The density-heating rate relationship is valid for a single loop. For a multi-thread 
simulation we assume the heating rate for each thread is related to this peak heating rate by 



where an determines the duration of the heating envelope. The parameter to is chosen so 
that all of the times are positive. The heating events are spaced so that as the heating in the 
previous loop ends the heating in the next loop begins. This heating scenario is illustrated 
in Figure [13] for aH = 100, 200, and 300 s. 

Once the individual hydrodynamic simulations are run, we average over the loop apex 
at each time step to compute a representative temperature and density. These densities 
and temperatures are then used as inputs to the TRACE temperature response to calculate 
the expected count rates in the TRACE 171, 195, and 284 A channels as function of time. 
The simulation times are shifted so that the peak in the 171 A emission corresponds to the 
observed peak. Since we are not interested in resolving the diflFerences in absolute calibration 



Ei{t) = Eo + git)Ei 



(9) 



{tj - to) 
2a% 



21 




exp 



(10) 
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among the various instruments, we also introduce a scaling factor so that the peak simulated 
emission matches what is observed. The simulation results are shown in Figure [HI and 
indicate that gh = 200 s simulation, which yields a simulated loop lifetime of at = 286 s, 
best matches the observations. 

We also use the simulated densities and temperatures as a function of time to compute 
the expected intensities in many of the emission lines that can be observed with EIS. Light 
curves from selected emission lines are shown in Figure [151 Since the absolute time for 
the simulation has been established through the comparisons with TRACE ^ we select the 
simulated EIS intensities that correspond to the time of the EIS observations and use them 
as inputs to the same differential emission measure code that was used to produce Figure [5l 
The resulting simulated DEM is shown in Figure [151 The simulated intensities are given in 
Table [H 

The agreement between the observed and simulated differential emission measure is 
relatively good. The simulation captures the salient features of the observations, a relatively 
high density and a narrow DEM. For these simulated intensities we do obtain a somewhat 
lower electron density (log rig = 9.5 for the simulation and 9.7 for the observation) and peak 
temperature (logTg = 6.04 for the simulation and 6.11 for the observation). The dispersion 
in the DEM also does not match the observation exactly (loga^ = 5.24 and 5.48). Given 
the approximate nature of the simulations we consider these discrepancies to be small. The 
difference in the density comes about because we have inferred the heating rate from a family 
of single-loop hydrodynamic simulations but the emission is actually a composite from several 
threads. It is likely that iterating on this solution would yield better agreement with the 
observations, but this is unlikely to yield addition physical insights. 

Finally, we have simulated the expected XRT emission for the Open/Ti-Poly filter com- 
bination using the standard XRT software routine xrt_t_resp. The simulated and observed 
light curves are shown in Figure [161 This comparison presents the greatest challenge to the 
modeling. The modeled composite intensities, which have been scaled to match the observa- 
tions at 20:00 UT, clearly do not extend back in time enough to cover the entire evolution of 
the emission in this region. The peak observed emission in this loop occurs at approximately 
19:30 UT, before the simulation has even begun. 

Given the diffuse nature of the XRT emission and the difficulty of isolating individual 
loops at high temperatures perhaps the simplest explanation may be that the XRT light 
curve includes the contributions of many loops in addition to the loop we isolated using 
the EIS and TRACE data. The identification of individual loops at very high temperatures 
with XRT is likely to be hampered by the slow evolution of plasma during the conductive 
cooling phase. This is evident by the slow evolution of the threads in the simulation. In 
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Figures [15] and [161 for example, we see that the threads last for approximately 1 hour at 
high temperatures. At the lower temperatures, when radiative losses are much higher, the 
cooling is dominated by radiation and the evolution is much faster. In these simulations 
the threads last only for about 10 minutes in the TRACE 171 A bandpass. This difference 
suggests that the differentiated loops seen at lower temperatures, such as those illustrated 
in Figure [6l would appear as a single structure in XRT. These differences are also related to 
the broad temperature response of XRT. The relatively narrow line emission imaged with 
EIS and TRACE emphasizes small differences in temperature. 

Alternatively, the inability of the model to reproduce the observed XRT emission may 
reflect inadequacies with the hydrodynamic simulations during the conductive phase of the 
cooling. It may be that the heating is not as impulsive as we have assumed. Many previous 
studies have suggested a gradual decay in the heating (e.g. jReale et al.ll2004l ). These differ- 
ences in the assumed heating may be related to the changes in the topology of the magnetic 
field during the evolution of the event. Observations of post-flare loop arcades have shown 
that newly reconnected field lines re lax from cusp-shaped to approximately semi-circular dur 



ing the early phases of t he cooling ( jSvestka et al.lll987l : iForbes fc Actonlll996l : ISheeley et al 



20041 : iReeves et al.ll2008l ). The comparisons between the hot and cool emission shown in Fig- 
ure [10] clearly suggests that field line shrinkage is occurring in this event. The cool post-flare 
loops observed with TRACE are generally observed at the lowest heights of the arcade and 
do not overlap with the high temperature XRT emission seen at the top of the arcade. The 
implications of field line shrinkage on hydrodynamic simulations has not been investigated. 
In general, the conversion of magnetic energy into thermal e nergy through the pro cess of 
magnetic reconnection is not well understood (see, for example, iLongcope et al.ll2009l ). More 
detailed analysis of MHD sim ulations is needed to better understand the evolution of coronal 
loops after reconnection (e.g.. iLinton fc Longcopdl2006l ). 



5. Summary and Discussion 

In this paper we have made use of the unprecedented opportunity to observe evolving 
coronal loops in detail. We have used STEREO^ EIS, TRACE^ and XRT data to constrain 
a multithread model of coronal heating and compare with observations. These comparisons 
indicate that it is possible to reproduce the high densities, long lifetimes, and relatively 
narrow emission measure distributions inferred from the data so long as the heating envelope 
of the heating is sufficiently narrow. 

The most challenging comparisons are with XRT, where the model fails to reproduce 
the extended lifetime of the emission at high temperatures. It is not clear if this is due to 
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our inability to isolate narrow loops at high temperatures or to problems with the assumed 
envelope on the heating. Recent analysis of the EIS spectral range has identified Ca XIV, 
XV, XVI , and xvil emission lines that can be used in the analy sis of high temperature 



plasma (jWarren et al.ll2008l : iDel Zannall2008l : IWatanabe et al.ll2007l ) . These lines will provide 



additional information on plasma evolution during the conductive phase and have been 
incorporated into the latest EIS observing sequences. New active region observations should 
be available during the rise of the next solar cycle. 

Ultimately our goal is to apply the multithread modeling described here to non-fiaring 
active region loops. It is encouraging that the simulated EIS diflPerential emission mea- 
sure curve derived here is similar to those derived for lower density active region loops 



( IWarren et al.ll2008l ). It remains to be seen, however, that this models can also reproduce 
the loop evolution observed in TRACE and XRT. The launch of the Atmospheric Imaging 
Assembly (AIA) on the Solar Dynamics Observatory (SDO)^ which will combine full disk 
imaging, TRACE-like spatial resolution, 10s cadences, and multiple filters, will greatly ex- 
pand the number of useful active region observations that combine EIS plasma diagnostics 
and loop evolution. 
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Fig. 1.^ A summary of the GOES, EIS, TRACE, XRT, and EUVI data taken on May 2, 
2007 between 17:00 and 23:00 UT. Each image is indicated with a dot. For EUVI the times 
of the STEREO A images are shown below the STEREO B images. There are gaps in the 
TRACE data due to orbital eclipses. There is also a gap in the STEREO A data. 
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Fig. 2. — EIS rasters of active region 10953 in 10 different emission lines. These images were 
constructed by stepping the V' sht from west (right) to east (left) over the region between 
20:27 and 21:38 on May 2, 2007. The field of view is 256'' x 256''. Each exposure is 15 s in 
duration. Note that the Fe XI 192.813 A raster shows both Fe XI and Ca XVII 192.858 A 
emission. The dark vertical band is from a Hinode orbital eclipse. For comparison with 
the imaging data a vertical line has been drawn on the rasters corresponding to data taken 
20:58:34 UT. 




Fig. 3.— EVYl/ STEREO B 284, 195, and 171 A images of active region 10953 from ap- 
proximately 20 to 22 UT on May 2, 2007. The field of view is approximately 216" x 168". 
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Fig. 4. — Selected XRT Open/Ti-Poly images from active region 10953 from 17 to 23 UT 
on May 2, 2007. The field of view is 281" x 226". 
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Fig. 5. — An EIS snapshot of an evolving loop. The top left panel shows the region surrounding 
the loop in Fe XII 195.119 A. The pluses (+) are the spline knots used to define the coordinate 
system along the loop. The bottom panels show the straightened loop region in s^t coordinates 
for 9 emission lines. The intensity averaged along the loop is also plotted for each line. The 
region used to compute the background subtracted loop intensities is indicated with the plus (+) 
symbols. For reference the Fe XII 195.119 A intensities are also shown in each panel (dotted line) 
and the Fe XII 195.119 A centroid is also shown in each panel (vertical line). Finally, the correlation 
between the Fe XII 195.119 A intensities and the loop intensities for each line are indicated. The 
loop emission is well correlated at temperatures below Fe XIII and poorly correlated at higher 
temperatures (Fe XIV-xvl). This temperature dependence is reflected in the differential emission 
measure shown in the top panel. Here the emission measure loci are shown (black lines) along with 
the DEM (red line). 
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Fig. 6.^ EIS multicolor images composed of Fe XVI 262.984A (red), Fe XII 195.119A 
(green), and Si VII 275.354 A (blue). Each wavelength is used as one of the color channels in 
the composite 24 bit image. The panel on the left shows the active region, the panels on the 
right show the loop. The panel on the far right is the composite loop image. In regions where 
there is strong emission at multiple temperatures the image is white. The loop arcade is 
dominated by the primary colors suggested a relatively narrow distribution of temperatures 
the loops. 




Fig. 7.— TRACE 171 A intensities as a function of time for the loop observed with EIS. 
The TRACE and EIS images near the peak of the event have been co-ahgned. The vertical 
lines indicate the times for the TRACE images shown in the figure. 
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Fig. 8. — STEREO EUVI reconstruction of the loop geometry. Because of a data gap with 
STEREO B, images taken May 2, 2007 21:16 UT, shghtly after the loop appears in TRACE, 
are used. The top panels show the loop traced out in both images. The image on the far 
left is the STEREO B 171 A image taken near the peak in the TRACE 171 A light curve. 
The bottom panels show the loop in various planes. The loop length is 135 Mm and the loop 
inclination is 68.5°. 
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Fig. 9. — XRT Ti Poly intensities as a function of time for the region observed with EIS. 
The vertical lines indicate the times for the XRT images shown in the figure. The crosses 
indicate positions along the loop used to define the loop coordinate system. 
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Fig. 10.^ Multicolor images from XRT and TRACE. The red channel is TRACE and the 
green channel is XRT. The image times are offset by one hour to illustrate the differences in 
morphology of post-flare loops in different stages of evolution. The high temperature XRT 
images do not show the fine loops observed at lower temperatures in TRACE. 
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Fig. 11. — Left panels: XRT and TRACE images of the loop at different times and temper- 
atures. The region used to compute the loop intensities are indicated. Right panels: The 
average XRT and TRACE intensities across the loop. The dotted line indicates the assumed 
background. 
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Fig. 12. — A family of hydro dynamic simulations for the loop geometry inferred from 
STEREO. The top panels show the apex densities and temperatures as a function of time. 
For each simulation we have determined the density as the loop cools through 1.3 MK. The 
bottom panel shows the relationship between the input heating rate and the density at 
1.3MK. 
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Fig. 13. — The heating rate for multi-thread simulations of the loop. The Gaussian curve 
represents the heating rate envelope. The triangles represent the heating rate for each thread. 
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Fig. 14.^ Simulated and observed TRACE light curves for an = 100, 200, and 300 s. 
The simulation times have been adjusted so that the peak emission in 171 A matches what 
is observed. The magnitude of the simulated emission has also been scaled to match the 
observations. The lifetime of the loop is calculated from a Gaussian fit to the composite 
light curve. The an = 200 s simulation most closely matches to observed loop lifetime. 
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Fig. 15. — {top panel) The apex temperature as a function of time for each thread. The 
vertical hne represents the approximate time of the EIS observations, {middle panel) Simu- 
lated EIS light curves for selected emission lines. Only the composite light curves are shown. 
{bottom panel) The DEM inferred from the simulated EIS emission. 
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Fig. 16. — Simulated and observed light curves for the XRT Open/Ti-poly filter combination. 



Table 1. Observed and Modeled 
Intensities^ 



Line 


^ohs 


-^dem 




-^dem 


Si VII 275.352 


37.2 


37.4 


66.7 


61.2 


Fe X 184.536 


190.3 


194.4 


339.5 


342.6 


FeXI 188.216 


436.9 


357.9 


439.7 


496.7 


FeXII 195.119 


718.6 


853.2 


718.6 


773.7 


Fe XIII 202.044 


595.5 


545.1 


301.2 


277.1 


Fe XIII 203.826 


207.9 


200.3 


120.7 


120.9 


Fe XIV 274.203 


0.0 


106.1 


33.9 


28.7 


FeXV 284.160 


0.0 


296.1 


36.9 


31.2 


Fe XVI 262.984 


0.0 


2.71 


0.1 


0.1 



^Units are erg cm ^ s sr ^. Isim refers to 
the simulated intensities presented in Section [H 



